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Abstract 

We give a new elementary proof of the fact that the value of the least fc*'* power non- 
residue in an arithmetic progression {bn + c}„=o,i...j over a prime field Fp, is bounded by 
7/-\/5 • b ■ \fpfk, + 4& + c. Our proof is inspired by the so called Stepanov method, which 
involves bounding the size of the solution set of a system of equations by constructing a non- 
zero low degree auxiliary polynomial that vanishes with high multiplicity on the solution 
set. The proof uses basic algebra and number theory along with a determinant identity that 
generalizes both the Sylvester and the Vandermonde determinant. 

1 Introduction 

Let Fp be the prime field with p elements. An element a G Fp is called a k^^ power non-residue 
if there is no 6 G Fp such that = a. Bounding the value of the least fc*'' power non-residue in 
a prime field Fp, where k \ p — 1, is a fundamental problem in number theory and algebra. It 
has an important application in finding roots over finite fields. For instance, it is known from 
the work of Vinogradov [Vin72| (see also Proposition 7 in |Evd94| ) that given a k^^ non-residue, 
all the k^^ power roots of an element a G Fp i.e. all x such that = a, can be found in 
(k ■ logp)*^^^^ time. It is a major open problem in number theory to show that the least A;*'^ 
power non-residue is bounded by {logp)^^^\ Indeed, such a bound is already known under the 
powerful assumption of the Extended Riemann Hypothesis (ERH). It follows from the work of 
Ankeny [Ank52] and Bach |Bac82j that assuming ERH , the value of the least A;*'' non-residue 
in Fp is bounded by 0(log^ p), where /c is a prime dividing p — 1. However, such a strong bound 
is not yet shown without the assumption of any unproven conjecture. We now briefiy mention 
the known results on ERH-free bounds on least power non-residues. 

1.1 Earlier work 

The Polya- Vinogradov inequality (see Chapter 23 in [DavOOj ) states that 



E 



m+l<x<m+n 



< VP - logp. 



where x is a non-principal character modulo p. Taking x to be the quadratic character, it 
immediately follows that the least quadratic non-residue in Fp is bounded by y/plogp. In 
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1919, this bound was improved by Vinogradov (see [Vin54irVin85j ) . who showed that the least 

quadratic non-residue in ¥p is less than p^V^ \og^ p. In a subsequent work, Vinogradov |Vin27] 
also showed that if k \ p — 1 and k > m^, where m is an integer greater than 8, then the 
least A;*^ power non-residue is less than p^/™ for all sufficiently large values of p. Later, in 
1957, Burgess |Bur57j improved upon Vinogradov's result and showed that the least quadratic 

non-residue is in fact bounded by p'^V^'^'^ for any small enough e > 0. A simple account of 
Burgess' theorem can be found in the work of Stepanov |Ste75| (see also [Kar68] ). We note 
that the proofs of Vinogradov and Burgess' results involve sophisticated analytic arguments on 
character sums. On the other hand, using purely elementary methods Brauer [Bra32j showed 
that the length of the largest sequence of consecutive A;*'* power residues or non-residues is 
bounded by ^/2p + 2. Later, Hudson [Hud74| gave an elementary argument to show that 
the value of the smallest k^^ power non-residue in an arithmetic progression {bn + c}„=o,i... is 
bounded by 211/455/2^2/5 ^g^3y/5 ^26^ if p is sufficiently large. Surely, these bounds are worse 
than the best known bounds of Burgess and Vinogradov. Nevertheless, it is perhaps interesting 
to know how much elementary methods can achieve in proving non-trivial bounds for power 
residues and non-residues. 

1.2 Our results 

We give a simple proof of the following fact. 

Theorem 1.1. The value of the least k^^ power non-residue in an arithmetic progression {bn + 
c}jj=o,i,... over Fp is hounded by 7/ y/b ■ b ■ ^ + Ab + c. 

Notice that, for k > p^^^, the bound given by Theorem II. II is better than the bounds shown by 
Hudson |Hud74] and Brauer |Bra32j . Our proof is inspired by the polynomial method, which 
was introduced by Stepanov to give elementary proofs of many of the significant special cases 
of Weil's theorem on rational points on curves. The reader is encouraged to refer to the book 
by Schmidt |Sch04j for an account of the elementary methods used in studying equations over 
finite fields. (For a quick introduction to some of the main results in this area refer to Tao's 
blog entry [EqQ9].) 

The main idea behind Stepanov's method is to construct a non-zero auxiliary polynomial 
that vanishes with high multiplicity on the solution set of a system of equations. Now, if the 
degree of the auxiliary polynomial is also 'small' then this can be used to upper bound the 
size of the solution set. We use this theme of bounding a solution set size via a low-degree 
auxiliary polynomial to give a new elementary proof of the square root bound on the least 
power non-residue/residue in any arithmetic progression. But, it turns out that the only 'not 
so easy' part of our proof is showing that the auxiliary polynomial thus constructed is non- 
zero. We resolve this difficulty by using an interesting determinant identity that generalizes 
the determinant of both the Sylvester and the Vandermonde matrix. Proving this determinant 
identity constitutes the main technical contribution of our work. We hope that this identity on a 
generalized Sylvester- Vandermonde matrix is of independent interest and may find applications 
elsewhere. 

2 The Polynomial Method 

In this section, we describe our approach to proving Theorem ll.il At the heart of our argument 
is the following lemma. 
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Lemma 2.1. A system of univariate polynomials S = {{x + Oj)* — 0t}i<j<r, where Oi,ai G Fp 
and Qi 's are distinct, has at most 2t/{r — 1) +3 common roots, if r < 2/\/h- y/t + 1 and p > 2t. 

Before we prove this lemma, let us at first see how it implies Theorem 11.11 (To keep the 
presentation simple, we avoid the use of the floor/ceiling notations. The analysis can be made 
more precise, at the cost of making the constant 4 in Theorem II . II and the constant 3 in Lemma 
12.11 slightlv worse.) 

In Lemma 12.11 take Oi = b ■ {i — 1), t = {p — ^)/k and 9i = 1 for all 1 < i < r. Set 
r = 2/^^h■^^t + l. If the sequence of elements hj + c,bj + c + h,hj + c+2h,...,hj + c+h{r — l) are 
A;*^ power residues then surely, hj+c is a common root of the system S. By Lemma [2.H there are 
at most 2t/(r — 1) + 3 common roots of S. In the worst case, all these common roots can possibly 
be consecutive elements of the arithmetic progression {bn + c}n=o,i,...- Therefore, the first index 
m for which bm + c is a k power non-residue, can be at most 2t/(r-l) + 3 + r = 7/^fh-y/t + A. 
The same argument can be used to prove a slightly general form of Theorem 11.11 as stated in 
the following corollary. 

Corollary 2.2. The length of the largest sequence of consecutive k^^ power residues or non- 
residues in an arithmetic progression {bn + c}„=o,i,... is bounded by 7/-v/5 • ^Jp—\lk + 4. 

The rest of this section and the following section (Section [3]) are devoted to the proof of Lemma 
2.1 Proof of Lemma 12.11 

The strategy we employ to bound the number of common roots of 5, denoted by z^(5) henceforth, 
is inspired by what is known as the 'Stepanov method' (also called the 'polynomial method'). 
The idea is to show the existence of a non-zero polynomial F{x) of small degree (say) such 
that if a is a common root of S then a is also a root of F(x) with multiplicity M (say). If this 
happens then we immediately know that i>{S^ can be at most By making N as small as 
possible and M as large as possible, we can arrive at an upper bound for z^(5). 
Let us see how to put this idea at work. Choose F{x) to be of the form, 

r 

F(x) = 5^G,(x)(x + a,)*+^^-l+^ 

i=\ 

where Gj's are polynomials of degree at most d (say) and M is the multiplicity parameter 
mentioned above. The parameters d and M will be fixed eventually in terms of t and r. Define 
s as, 

s = 0, if (r-1) I (t + M-1) 

= (r — 1) — ((t + M — 1) mod (r — 1)), otherwise. 

The role of the parameter s is to make t + M — 1 + s perfectly divisible by r — 1, a technical 
requirement for the analysis in Section [3] to go through. Let us take a short digression and 
clarify a bit more the purpose of the parameter s. 

One might wonder as to why we do not assume, for the sake of simplicity, that s = and 
t -\- M — 1 is divisible by r — 1. At some point in our argument we need to establish linear 
independence of a certain linear system. If the coefficient matrix associated to the linear system 
is a square matrix then all we need to show is that the corresponding determinant is non-zero. 
It turns out, it can be shown that such a determinant is non-zero by using an identity involving 
derivatives of the determinant function. Whereas, for a non-square system it is a little more 
tedious. 
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Coming back to the main flow of the proof, let us see what is required from the polynomial 
F{x). Denote the i^^ derivative of F with respect to x by F^^'> and let T = t + M — 1 + s. Also, 
G^^^ denotes the j*'^ derivative of G. If a is a root of the system S then we require F^^^ (a) = 
for all < £ < M — 1 since we want a to be a root of F with multiplicity M. This means, 

F^'Hc^) = E E ^^•(^) • • ^?^(«) • (" + ad''^'^'-^'"'^ = 0, (2) 

i=l j=0 

where Cj{T) = Y[t^o~^ — A;) is a constant and (a + Oj)* is evaluated to 9i since a is a root 
of S. Suppose that the coefficients of the polynomials Gj's, in Equation [U are variables. Also, 
treat the expression given in Equation [2] as a polynomial in a of degree {d + M— 1 + s — i) 
with coefficients as linear forms in the variables (that are coefficients of the Gi's). By equating 
these coefficients to zeroes, we can ensure that FW(a) is zero. Therefore, for any particular 
i, Equation [5] imposes {d + M + s — i) homogeneous linear constraints, yielding a total of 
M ■ {d + s)+ ^^(^^+^) homogeneous equations in (d+l) -r variables (as £ runs from to M — 1). 
Thus, in order that we get a nontrivial solution for the coefficients of Gj's, it is sufficient to 
satisfy the the following condition, 

M-{d + s) + ^^^^^^<{d+l)-r. (3) 

Further, we also need to ensure that this solution is such that F{x) ^ 0. The degree of the 
polynomial F(x) isN = d + t + M — 1 + s. If the number of variables (d + 1) • r is greater 
than d + t + M + s then surely there is a nontrivial setting of the coefficients of Gi 's that makes 
F{x) = 0. However, such a situation can be possibly averted if we also put the restriction that 

{d+l)-r<d + t + M + s. (4) 

Indeed, we show (in Section [3]) that Condition H] is sufficient to guarantee F{x) 7^ if the 
coefficients of the Gj's are not all zeroes. To summarize. Condition [3] ensures that we are able 
to find nontrivial Gi 's by solving the homogeneous linear equations arising from Equation [2l 
for < i < M — 1. Whereas, Condition [5] guarantees that the polynomial F{x), defined in 
Equation [H is non-zero if not all the Gj's are zeroes - the proof of this appears in Section [3l 

Putting together Condition [3] and [H and using D = d+l, we get the following overall 
condition to satisfy. 

M • (Z? + s) + <D-r<D + t + M + s-l. 

Since our objective is to minimize the quantity ^ = ^ would like to minimize 

D + t + M + s — 1, which being lower bounded by D • r, the best we could possibly do is to 
choose D such that, 

D-r = D + t + M + s -1 

= i±^i±l^ (5, 

r — 1 

This setting of D satisfies Condition HI Now, let us see how to satisfy Condition [3l Choose 
M = r/2 and put D = (t + M + s — 1) /(r — 1) as in Equation [5j Using the fact that s < r — 1 
and then simplifying further, Condition [3] reduces to the following quadratic inequality: 
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It is easy to check that this is satisfied if r < ■ \ft + 1. We are almost done. Recall 

that the maximum size of v{S\ the set of common solutions of 5, is bounded by N/M , where 
N = deg(F). Hence, 

N D + t + M + s-l D-r , ^ 
W{S)\ <j^< = (usmg Equation 15]) 

Since M = r/2, \i'{S)\ < 2D. Once again, using the value of D it is easy to derive that 

IK5)|<^ + 3. 
r — 1 

This proves Lemma l2 . 1 1 except the lemma: 

Claim 2.3. If D = ^^y+j'"^ then F{x) = if and only if Gi{x) = 0, for alll<i<r. 

The next section is devoted to the proof of this statement. The main ingredient of the proof is 
an identity involving a generalized Sylvester- Vandermonde determinant. The condition "p > 2t" 
(in Lemma l2.ip also appears in this proof. 



3 A Generalized Sylvester- Vandermonde Determinant 

Recall, from Equation [H that F{x) is defined as F{x) = Y^\=i^i{^) ' + ai)^ , where T = 
t + M + s — 1. Suppose Gi = l^j=o CijxK Then, 

r d 

= ^^^^Cijx\x + aif. 

Proving Claim 12.31 essentially means proving this: if F{x) = then Cij = 0, for all i and 
j. Suppose, on the contrary, that this is false. Then, the polynomials {x^{x + aj)-^}jj', for 
1 < i < r and < j < d, are F-linearly dependent. In other words, the following matrix. 





■■ 


• i/^dK-' ••• 


••• CoK 








CM ••• 


••• CM •• 


• cx 






\T-dl "2 


••• m 




V = 




Cr)al ■■■ 


••• CK •• 


■ CK 










. 1 T \ T-d ... 

\T-dJ 


••• CK 






V 




■ CK •• 


■ CoKJ 



must be singular. But, we show, in Lemma l3.ll that V cannot be singular if a^'s are distinct 
and p > 2t. This leads us to the necessary contradiction and hence a proof of Claim (231 



Remark - Notice that, V is a square matrix since the number of rows D ■ r equals the number 
of columns D + T, by the choice of D = T/(r — 1) in Claim 12.31 We call V a generalized 
Sylvester- Vandermonde matrix because when r = 2, it becomes the Sylvester matrix of the 
two polynomials {x + oi)^ and {x + 02)"^, whereas when D = 1, it is the Vandermonde matrix 
(scaled appropriately) . 
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We now prove the following identity. 
Lemma 3.1 (Sy I vaster- Vandermonde identity). The det(V) = C ■ Y[i<i<j<r ~ '^j)^^ > where 

It is not hard to check that C 7^ in Fp, if p > 2t (just use the facts that s < r — 2 and r is an 
integer less than or equal to 2/-v/5 ■ \/t + \). 

3.1 Proof of Lemma 13.11 

First, we show that det(V), viewed as a polynomial in ai, . . . , a.^, is divisible by (ai — a2)^^ ■ 
Then, by symmetry, det(V) is also divisible by {ai — aj)^ for every pair (i, j) with i < j. Hence, 
det(V) is divisible by ni<i<j<r ~ ^j)^ • looking at the matrix V, it is easy to infer that 
the highest degree of oi in det(V) (once again, viewed as a polynomial in ai, . . . , ar) is at most 
D ■ T. Since the degree of ai in the expression ni<t<j<r ~ ^i)^ ■ {r — I) = D ■ T {as 

D = T/(r — 1)), det(V) must be of the form C -YiKi^j^r i^^i ~ ^3)^^ ■> where C is just a function 
of r, d and r, but not the Oj's. 



In the proof, it will be more convenient if we express matrix V in terms of polynomials. 
Notice that the rows of V can be identified with the coefficient vectors of the polynomials 



{x + aif , x{x + aiY , . . . ,x'^{x + ai)'^, (j; + 02)^, x(x + 02)^, 
us abuse notations slightly and write V as, 

x{x + ai)"^ 



V 



{x + a2Y' 
x'^{x + 02)-^ 



{x + arY' 



\x"-{x + ttr) I 



. ,x'^{x + a2) 

( (x + ai)^ \ 
(x + ai)^+i 



(x + ai)'^+'^ 



and so on. Let 



row operations 




(x + arY 



\{x + a,)^+V 



(6) 



meaning that V is formed by the coefficient vectors of these polynomials. Let Rjj be the row of 
V standing for the coefficient vector of x^{x + aiY' . Consider the following row operations on 
V. 

Ry ^ X] (i) ■ ■ ^^i-^^ for 1 < i < r and < i < 
fc=o ^ ' 

Equivalently, after the row operations, the coefficient vector of x^ix^aiY' gets replaced by that 
of the polynomial, 

Y.(^\ ■ a\ . x^'\x ^ = {x ^ a,f^K 

This leaves us with a transformed matrix V, as shown in Equation[6l such that det(V) = det(V'). 
To show that (01—02)^ divides det(V'), view det(V') = /(ai) as a polynomial in oi. Let f^^\ai) 
denote the order derivative of /(ai) with respect to ai. It is sufficient if we are able to show 
that ai = 02 is a root of f^^\ai), for all < £ < Z)^. 
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Claim 3.2. Let /(oi) = det(V') and /^(ai) = ^/(ai). Then 

0-1 =0,2 is a root of f^^\ai), 

for aUO<e< D"^. 



Proof. To prove this claim, we need the following identity involving derivatives of a determinant. 
Let A = {aij) be an n x n matrix whose entries are real functions of x. Then, 



\ d^" „ d^" „ d^" ^ 



where (^^ ^ ) multinomial coefficient. Now imagine applying this identity to f^^\ai), 

the ^^'^ order derivative of det(V'), where the entries of V are viewed as functions of ai. It is 
clear from Equation [6] that except for the first D rows of V, the rest are independent of the 
variable ai. Denote by Rjj, the row of V generated by the coefficients of {x + ai)"^~^^ . We write 
^Rjj to mean the row formed by applying the operator ^ to every entry of Rjj. Therefore, 
we have the following identity. (For economy of space, we switch to the transpose notation.) 

/ £ \ ( d^^ d^-o ^ 

/ («l) = 2^ [p P ) "i^* [tTT^IO, —^Rld, R20, R2d, RrO, Rrd] 

Notice one nice property of the polynomial representation of V: In the above equation, -^^-^Rij 
is exactly the row formed by the coefficients of -^-rqrfC^ + '^i)^'^"' = Cj • {x + a\)^^^~^^+^ ^ where 

da^ 

Cj is a constant (depending only on T, j and ij+i). Now let us see how large £ needs to be so 
that (ai — 02) does not divide f^^\ai). 

Suppose (ai — 02) f f^^\ai). Then there exist a term L in the above summation that is not 
divisible by (ai — 02). Let that term be identified by some tuple {£i,£2, ■ ■ ■ ,£d)- Observe that 
this term, 

L = det [— j-Rio, . . . , „ Rid, R20, ■ ■ ■ , R2di ■ ■ ■ , RrO, • • • , Rrd] 

y da^ da\° j 

= det {[co{x + aif-'\. . . , Cdix + aif+''-'^,{x + 02)^, . . . , (x + 03)^+^ . . .] 

If any of the exponents {T — £i,T + 1 — £2, ■ ■ ■ ,T + d — £d} is greater or equal to T then 
(ai — 02) I L; since otherwise some row Cj{x + ai)^^^~^^+^ becomes equal to some other row 
{x + a2Y'^^ (up to a multiple of Cj), when ai is replaced by 02- Also, if any two of the exponents 
{T — £i,T + 1 — £2, ■ ■ ■ , T + d — are the same then L = 0. This leaves us with only one option 
- the set {T - ^1, T + 1 - £2, • • • , T + d - £d} is 'dominated' by the set {T - 1, T - 2, . . . , T - 
(We say a set 5i is dominated by another set 5*2 if for every element ei G Si there is a unique 
element 62 G 5*2 such that ei < 62.) Therefore, 

d D 

- 2 2 

= (Taking d = D - 1) 

□ 
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It follows from Claim lHT^ and the discussion before that det(V) is of the form C''ni<i<j<r ~ ^i) 
where C is a constant depending only on T, d and r. What remains to be done, in order to 
complete the proof of Lemma 13.11 is to show that C = Yie=o C^'^) 1 11^=0 (^i"*^) • '^^^ proof of 
this is included in Appendix Rl 

4 Discussion 

Although, the square root bound of Corollary 12.21 is not the best known bound for this problem, 
it may be worthwhile exploring the 'polynomial method' further to see if the bound can be 
strengthened, or if nontrivial bounds of some other related problems can be derived through it. 
Towards this, we have the following three questions in mind. 

Question 4.1. (Strengthening Lemma [2l]) Is it possible to give a better bound on the number of 
common solutions of the system S (defined in Lemma [2.ip . perhaps by considering an auxiliary 
polynomial of the form. 



Note, in our case, the auxiliary polynomial F{x) is defined with £ = 1. 

Question 4.2. (Simultaneous quadratic character) Is it possible to use our approach to show 
that for any pair of distinct elements a,b £ ¥p, the value of the largest possible m for which 
x((a + i) • (6 + i)) = 1, for all < i < m, is 0{y/plogp)? Here x(o) denotes the quadratic 
character of a. 

Question 4.3. (Least primitive element in ¥p) Is it possible to show that the value of the least 
primitive element in Fp is ■ poly(logp)) using our approach? 
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A The constant in Lemma 13.11 

Since the monomial YVi=i has coefficient 1 in the product ni<j<j<r(^i ~ ^j)^^ (viewed 

as a polynomial in the fli's), the constant C in Lemma |3. II is the same as the coefficient of the 

monomial rii=i *^ det(V). 

Claim A.l. The coefficient of the monomial Y\l=i det(V) is Y\J^q C"^'^) 1 11^=0 ("^j"'^) 

Proof. By definition, det(V) = J^aeS^ sign(cr) -n^gj^] where m = Dr, Sm is the symmet- 

ric group of degree m, and Vij is the {i,jY'^ entry of V. Note that, every product H^efm] ^e.<j{e) is 
a monomial in the Oj's with an attached coefficient. We need to find out, which all permutations 

a give rise to the monomial n[=i 

Let both Ri and Ci denote the set {D{i -1) + 1, . . . , D{i - 1) + so that [Dr] = i?i U • • • U 
Rr = Ci L) ■ ■ ■ L) Cr- We think of the i?j's as partitioning the rows and the Cj's as partitioning 
the columns of V. For instance, the rows with indices in Ri contain only those terms involving 
the variable a^. A crucial observation here is the following. The monomial is 
generated by exactly those permutations a that induce bijections between Ri o Ci, for all 
1 < i < r. This gives us a strategy to find the coefficient of Y[i=i 

Define the matrix Mi as the D x D submatrix of V which is induced by the rows Ri and 
the columns Cj. Since each term in det(Mj) has the same degree in a^, which is D^{r — i), the 

coefficient of Yii o-f ^'^ *^ can be obtained from the product ni<i<r det(Mj). 
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Notice that det(Mj) = af • det(-ffj), where Hi is the fohowing matrix formed by the 
binomial coefficients of the terms in Mj. 



/ ( 



Hi 



T 

T-D(i-l), 

( ^ 
\T-D{i-1)+ 



VT-D(i-l) 



l) 



T 

T-D{i-1) 
T 



\ 



\{T-D{i-l)+d) (T-D(i-l)+(d-l)) 



\T-D{i-l)-{d-l) 



(T-D(i-1)) / 



Therefore, the coefficient of Y\ 



r D'^{r-i) 



in det(V) is exactly Y\^^-^j. det{Hi) . For 1 < i < r, 



each of the binomial coefficients in the matrix Hi is non-degenerate, whereas for i G it is 

easy to see that det{Hi) = 1 as Hi is a triangular matrix with units along the diagonal. Suppose 
1 < i < r. After an appropriate row transformation, Hi gets transformed to. 



H 



T+d \ 
T-D{i-l)+d) 

( T+d-1 \ 
\T-D(i-l)+d) 



T+d 



\T~D{i-l)+{d-l) 



1.) 



\T~D{i-l)) \ 

( T+d-1 \ 
\T-D{i-1)) 



\{T-D(i-l)+d) il 



\T-D{i- 



\T-D{i-~l)+dl \T-D[i-l)+(d-l)l 

SO that det(i/j) = det(i^j'). Now we apply the following lemma, which we prove shortly, to find 
an expression for det(ff-). 



Lemma A. 2. For n,m,i £ N, satisfying i + m < n, 

/n+m\ 

\ £ ) 

fn+m—l^ 



/n+m\ f n+m \ 

{i+mj \i+m-l) 

/n+m—l\ /n+m—l\ 

\ i+m ) \l+m-l) 



Taking n ^ T, m 



nm /n+m\ 
j=o [ e+j ) 

nm fn+m\ 
3=0 y j ) 

ie+m) ie+m-i) " ' ie) 
d and £ — )• T — D[i — 1), in the above lemma, we get 



(iei{Hi 



d 

n 

3=0 



( ^+'^ 1 
\T-D{i-l)+j) 



n 



( T+d 
W{i~l)+j 
/T+d\ 



(7) 



Note that the condition ^ + m < n, in Lemma IA.21 is satisfied after the substitution since 
1 <i <r. Also, the above formula [7] evaluates to 1 for i G {l,?*}- Hence, the formula holds for 
all 1 < i < r. Therefore, the coefficient of ni=i 



D2(r-i) . 

a.. IS, 



r r d ( T+d 

n<iet(ff,)=nn%Sf 

1=1 i=lj=0 \ j J 



r\T+d (T+d\ 
11^=0 



□ 



as claimed. 

It remains to prove Lemma |A.2[ 
Proof of Lemma \A.S[ Index the rows from the bottom - the O*'^ row Rq is the bottommost. 

\ I ) 



Consider the row operation (—1)™' • Rm + Yl 



m~l , 
k=0 ' 



the topmost row in the i*'^ column (from the right, starting from zero) is 



Ru. Then the value of 



E 

fc=0 



n + k 
£ + i 



m 



fn+m\ fn+m\ m 
\ £ ) _ \ e ) 



(£+i\ 

\ £ J k=0 



III / ^ 

1 — n V / 



n 



£ + k 
i 
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The claim is - the quantity inside the summation is zero for i S {0, 1, • • • , m— 1}. This is because 
of the fohowing identity involving binomial coefficients (see page-169 in the book |GKP89j ). 



which is zero as (j^^) = when i < m (by definition). For i = m, the above equation evaluates 
to (—1)™". This means, after the transformation the value of the leftmost entry of the top row 
is (— 1)™ • / J whereas all the remaining entries of the row are zeroes. Recall that, 

in the row transformation we have multiplied the first row R^, by (— 1)"*, and so remultiplying 
by (— 1)™- the top-left entry becomes simply ("^™) / (^^™) • Using this argument inductively on 
the minors, the determinant evaluates to, 

™ /n+m\ 

n \ I ) ^ TT V e+j ) 
(e-+j\ 1 i (n+m\ ■ 

j=o \ e J j=o V j ) 



The last equality is simple to verify. □ 
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